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The collective modes observed in the loop-current ordered state in under-doped cuprates by 
polarized neutron scattering require that the ground state is a linear combination in each unit-cell 
of the four basis states which are the possible classical magnetic moment configurations in each unit- 
cell. The direction of such moments is in the c-axis of the crystals. The basis states are connected by 
both time-reversal as well as spatial rotations about the center of the unit-cells. Several new features 
arise in the theory of polarized neutron scattering cross-section in this situation which appear not 
to have been encountered before. An important consequence of these is that a finite component 
transverse to the classical magnetic moment directions is detected in the experiments. We show 
that this transverse component is of purely quantum-mechanical origin and that its direction in 
the plane normal to the c-axis is not detectable, even in principle, in experiments, at least in the 
quantum-mechanical model we have adopted. We estimate the direction of the "tilt" in the moment, 
i.e. the ratio of the transverse component to the c-axis component, using parameters of the ground 
state obtained by fitting to the observed dispersion of the collective modes in the ordered state. We 
can obtain reasonable agreement with experiments but only by introducing a parameter for which 
only an approximate magnitude can be estimated. Approximate calculations of the form-factors are 
also provided. 



PACS numbers: 



I. INTRODUCTION 



Polarized elastic neutron scattering experiments^'^ and 
dichroic ARPES'^ have revealed that the pseudo-gap 
phase of the cuprates have a long-range magnetic order 
which breaks time-reversal symmetry without breaking 
translational symmetry of the lattice. It is a Q = stag- 
gered order with zero net-moment in each unit-cell. Its 
geometric arrangement is consistent with the order of a 
pair of oppositely directed fluxes due to current loops 
formed in the o-cu-o links in each unit-cell^. Classi- 
cally such an order has four possible domains as shown 
in Fig. [T] These domains are specified by the directions 
(±1, ±1), that the order parameter Jl makes with respect 
to the X and y-axes of the crystal. The order parameter 
is an anapol^'^ given by 



Ct 



dr(L(r) x r). 



(1) 



where L(r) is the magnetic moment in the unit-cell at 
the point r. Such an order, for any of the four possi- 
ble domains, has orbital magnetic moments L pointed 
in directions along or opposite the c-axis of the crystals. 
However, polarization analysis of the neutron scattering^ 
has shown that this is not true. The direction of the 
moments, interpreted according to the classic theory of 
polarized neutron scattering^"®, makes a large angle with 
respect to the c-axis^; the direction along the plane is not 
revealed due to the multi-domain nature of the crystals 
and/or the multi-domain nature of the order or as we 
will show here due to its quantum-mechanical nature. 

Polarized inelastic scattering has also discovere d^";^^ 
two branches of weakly dispersive collective modes in 



the same temperature region as the magnetic order and 
with an intensity as a function of temperature compatible 
with it. Such collective modes can only be understood 
as quantum-fluctuations of the observed order, just as is 
true for the one branch of collective modes in the trans- 
verse field Ising modeU^. In the paper preceding thisi^, 
hitherto referred to as I, we have introduced a quantum- 
mechanical model for the observed order and calculated 
the simplest quantum-mechanical ground state of the 
model as well as the collective modes. The ground state 
is a product over the unit-cells of a sum over the four 
classical configurations in each unit-cell, depicted in Fig. 
[T] and given with a particular choice of phases by 
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is a parameter which has been determined by fitting 
the calculated collective mode dispersions to the exper- 
iments. I ± 1,±1),; refer to the four configurations in a 
unit-cell i. As discussed in detail in I, and will be sum- 
marized below, one can make unitary transformations 
consistent with the symmetry of the problem, which in- 
troduce other operators in the ground state, and give 
a corresponding ground state wave-function, which in 
general is a linear combination of the basis states with 
complex coefficients. We show in this paper that po- 
larized neutron scattering from such a ground state re- 
quires a quantum-mechanical description of scattering of 
neutron of the quantum magnetic moments, whereas the 
traditional method considers the problem as a quantum- 
mechanical scattering of neutron from a classical mag- 
netic field due to the ordered magnetic moments (and in 



some cases their zero-point fluctuations which only cause 
Debye- Waller like corrections in amplitudes). We show 
that in this situation, the observed "tilting" of the mo- 
ments is a purely quantum phenomena. One can only de- 
duce by neutron scattering, or by any other experiment, 
the component of the magnetic moment along the c-axis 
of the crystal and that perpendicular to it but not the 
two orthogonal components of the latter. We compare 
various aspects of the experiment with our calculations. 
Aspects of the symmetry of the observations all appear to 
be well reproduced. The quantitative magnitude of the 
"tilt" and its variation with the Bragg- vector can only be 
reproduced by introducing a free parameter, which can 
only be estimated approximately. 



The basis states in ([2]) may be taken to be the eigen- 
states of the orbital magnetic moment operator L3 , with 
the 3-axis identified with the c-axis of the crystal. An 
important aspect of the problem is to specify the kinetic 
energy term which mixes the four basis states in a cell i 
to give ©. This is important because there is obviously 
no orbital moment vectors pointing in the plane at the 
four locations indicated in the figure ([Ij because current 
flow only in the planes is assumed. We will show that an 
operator with the right commutation rules for an angular 
momentum exists in the problem but its physical basis is 
a fluctuating current loop between the four oxygen atoms 
around a cu atom in each unit-cell. Such an operator 
occurs naturally in the microscopic theory of loop cur- 
rents in the cuprateSii^ and leads among other things to 
the marginal fermi-liquidi^ii^ in the normal state in the 
quantum-critical regime. 



This paper is organized as follows. In the next sec- 
tion, we summarize the quantum Ashkin- Teller model 
introduced in I and cast it in the basis of local angular 
momentum operators which are more useful to discuss 
neutron scattering. We will also review the transforma- 
tion properties of these operators to show that, given the 
quantum Ashkin- Teller model, the direction of the "tilt" 
in the x-y plane cannot be determined. In the following 
section, we discuss neutron scattering and its polarization 
dependence and show that besides the usual matrix el- 
ement for polarized neutron scattering, there exists also 
another matrix element due to the flnite extent of the 
current loops. In calculating the neutron scattering in- 
tensity, we flrst consider the moments as point objects 
at the four sites and subsequently improve the calcula- 
tion by considering the flnite extent of the current loops. 
This also allows us to estimate the form factors or the 
momentum dependence of the scattering at the Bragg 
vectors. 



II. LOOP-CURRENT MAGNETIC ORDER 

In I we have fully described the symmetries of the 
quantum Ashkin- Teller model with which the collective 
modes of the ordered loop current states are described 
and compared with experiments. To calculate neutron 
scattering, we will proceed in two steps. First, we will 
stay with the abstract representation given in I. To repro- 
duce only the correct symmetry of the magnetic order, 
the orbital moments have been represented as point ob- 
jects located in the centroid of the triangular loops of 
Fig. dl]). The locations are labeled by R^^a = R? + R-a, 
a = 1, ..4. Here i denote the lattice sites, R" is the posi- 
tion of the center of the unit cell and Ra is the relative 
position of each local moments in the unit cell. We take 
Ri = (r-o,ro), R2 = {-ro,ro), R3 = (-ro.-ro) and 
R4 — ('"Oi ~''o) with a value rp smaller than 1/2 the lat- 
tice constant. The 4 loop current states are labelled by 
the eigenvalues of az and Tz of the classical Ashkin- Teller 
(AT) model. At the classical level, for state |1, 1) there is 
a magnetic moment perpendicular to the copper oxygen 
plane pointing up located at Ri and another pointing 
down at R3. Also there are zero moments at R2 and 
R4. All the other 3 loop current states are can be ob- 
tained by sequentially rotating |1, 1) by 7r/2 and will be 
denoted by | - 1,1), | - 1,-1) and |1,-1). This label- 
ing is also consistent with the direction of the anapole 
vector J7. The four states (±1,±1) are of-course glob- 
ally orthogonal. We assume that this is also true of 
the four local states in any unit-cell (±l,±l)i, as also 
each of them between different cells. This is no different 
than, say, what is done with respect to the local moments 
formed from the collective degrees of freedom of fermions 
in itinerant anti-ferromagnets or ferromagnets for regions 
of frequency and momenta where their exceptions do not 
overlap much the incoherent fermion excitations. 

This is adequate to get the symmetries of the neutron 
scattering intensity for different momentum transfer Q 
and different measured initial and final neutron polariza- 
tions. But it is not adequate to give the form- factor and 
relative intensities at different Q and polarizations. Cal- 
culating such information is a very formidable task. In a 
second step, we will however attempt this in an approxi- 
mate way by introduce the wave- functions responsible for 
generating the orbital moments as well as the representa- 
tion of the kinetic energy in terms of current-operators. 
Besides quantitative verisimilitude, this affords some in- 
sight into the interesting new physics in the kinetic en- 
ergy terms. 

The Hamiltonian derived in I to calculate the collec- 
tive modes was written in the space basis of the states 
(±1,±1), and equivalently in the direct product basis 
with 4-dimensional vectors. The transformation from one 
to the other is. 



In our conclusions, we discuss also alternate ways of 
obtaining the "tilts" and show that they are not consis- 
tent with the qualitative features of the experiments. 
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(3) 



The classical Ashkin- Teller model in this basis is given 
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Figure 1: The four Possible "classical" domains of the loop ordered state are shown. In the classical ordered phase, one of these 
configurations is found in every unit-cell. 



by 



Hat^~}_^ [JiStS'^ + J^tTf + JiKf'Kf] (4) 



The quantum terms causing a transition between the four 
states which are rotated with respect to each other in the 
direction of fl by ±7r/2. are given in the same choice of 
gauge as (0) by 



iJQ = ^t(S'i+ri)+t'X" 



(5) 



where S'^,T'^,K'^^ are matrices in S'C/(4) space, specified 
in I, where their commutation rules are also given. The 
Hamiltonian Hat + Hq was used to derive the collective 
modes in I and to get the ground-state wave function [51 
The terms K}^ which causes rotation by ±7r, are unim- 
portant for our purposes here because they cause change 
in angular momentum by 2 and therefore do not couple 
to neutrons, which can change angular momentum only 
by 1 in the weak scattering limit. 

For calculating neutron scattering, it is more conve- 
nient to define a basis set given in terms of the orbital 
moment operators at the 4 location in a cell. Since the 
local moments are generated by orbital loop current they 
should be considered in the representation for spin 1. The 
three components of the effective moment will be denoted 
by L = {L^,L^, L^). Normally, one would represent L by 
the spin-1 representation of SU{2), i.e. by a three di- 
mensional representation, with eigenvalues say of ±1,0. 
But we have four states per unit-cell in the loop cur- 
rent model. As explained in the Appendix [B] the three 
dimensional representation is inadequate to the present 
case and a four dimensional representation of the spin 1 
states must be used. This representation is given by 
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(6) 



This representation is equivalent to a representation in 
the basis of 2 independent set of spin-1/2 operators. This 
is quite natural to use for the quantum Ashkin- Teller 
model. We will verify that they satisfy the SU{2) algebra 
[L*,L^] =i€'^^^L^. This is a reducible representation with 
two (orthogonal) L3 = states we denote as |0^) and 
|0^). The loop current basis states in each unit-cell in 
this representation have the following bases in each unit- 
cell: 



|1,1> = (|1>r,,|01)r, 

|-l,l) = (|Oi)R,,h 
|-1,-1) = (|-1)r, 

11,-1) 



(|0'>R,,|1>R: 



|-1>R3,-|0^>rJ (7) 

1)r„-|02)r3,|1)rJ (8) 
|02)r.,|1)r3,-|0^)rJ (9) 

„-|01)r3,I-i)r4) (10) 



Here \i)n.j stands for the eigenvector of L3 with eigen- 
value i = ±1 or at location j = 1, ...4 in a unit-cell. The 
phase factors ±1 in front of the states |0)r. are picked 
so that the states in the left of Eq. ([7]) have zero net 
moment in a cell. 

While, since it has a classical analog, it is perfectly 
clear what Lf ^ physically means, the same cannot be said 
of L\ ^ and Lf ^ . Obviously, these are not proportional to 
angular momentum operators at the sites (i, s), since the 
currents are required to fiow only in the plane. We shall 
show that there exist operators which have off-diagonal 
matrix elements in the basis (±1,±1) , so that using 
([7]), we shall find (in general complex) matrix elements 
between the eigenstates with eignevalues ±1 and of the 
operator Lf^. We shall define L],. and Lf ^ through such 
matrix elements. 



A. Allowed Unitary transformations on the 
Quantum Ashkin-Teller Model 

The unitary transformations described in I are equiv- 
alent to rotations in L-space, as we will explicitly show 
below. Apart from the crystalline z-axis which defines 
the direction of the moments of the classical problem, 
there is nothing in the Hamiltonian derived in I for the 
problem which involves the crystalline x- and j/-axis (ex- 
cept that they be orthogonal to the z-axis). Therefore 
the direction 3 of L^ may be identified as the z-axis leav- 
ing the choice of the 1 and 2 axis with respect to the 
X and y crystalline-axes undetermined. This situation 
is similar to the traverse field Ising model in which the 
transverse field is self-generated. One can then choose 
it any direction normal to the Ising axis with identical 
experimental results. (Only if the transverse field is an 
externally applied field can one find in experiments that 
the direction of the moments is tilted from the Ising axis 
towards the external field axis.) 

We have shown in I that the classical Ashkin-Teller 
model has a continuous symmetry C/(l)s3 x U{\)t^ x 
[7(1)^33 which is the rotation around operators S^ , T^ 
and K^^. The rotation matrix is given by 



U{0) 



jeiS'y2ie2T^/2ie3K'^'y2 



(11) 



The classical AT model is invariant under this transfor- 
mation. It has been shown that the quantum terms in- 
troduced depend on the transformation U. At the same 
time, this transformation changes the wave functions. 
For a general wave function \ip) = a\l, 1) +b\l, —1) +c| — 
1, 1) -|- d| — 1, — 1) with four complex coefficients a, b, c, d, 
we have 



\U^) 



J4>3 



'D 



+ce"t"\~ 1,1) +de"^''\- 1,-1) (12) 



with 



1-1-92 



and 



^^ _ -61-02+03 ^ Thus U puts different phase factors on 
each of the 4 classical states with the sum of all phase 
factors restricted to zero. 

One should note the special case 9i = 62 = 0, when 
this unitary transformation is simply a rotation of S and 
equivalently of L in the x-y plane: 



UL^U'' = cosOsL^ - .-,111^3^ 
t/L^C/t = singgL^ +COS6I3L 



sin6^3L^ 
2 



(13) 
(14) 



So, given the quantum Ashkin-Teller model, the direc- 
tion of the moments perpendicular to the z-axis is un- 
determinable. 

We now show more generally that L^ have the prop- 
erties of angular momentum operators. As already dis- 
cussed ini^, by making use of Eq.([7]) etc, one can express 
the local spin operator L in terms of the Ashkin-Teller 



model operators of 5*' , T* (defined in^-) as follows 

-^R2 ^ -^R4 = 2'-'^ ^ ^ )' 

iR2=-iR4 = ^(^'-^') 

Using the commutation relations given for S'',T^ in I, it 
is easy to verify that the L satisfy the SU{2) algebra (up 
to an overall minus sign). Therefore, L can be regard as 
an angular momentum operators. 



III. POLARIZED NEUTRON SCATTERING 



The neutron scattering Hamiltonian is 



Hi, 



drB(r) • o-(r). 



(15) 



cr(r) is the spin of the neutron at point (r) in the crystal 
where the magnetic field operator is B(r). The source 
of the magnetic field are the magnetic moment operators 
due to spin or in our case orbital moments L(R^ + Rq) 
at locations (R^ -I- Ra). One can Fourier transform ([15]) 
and rewrite it in terms of the magnetic moments L at the 
momentum transfer Q as^'* 



^e'Q-^"F,(Q)L,a(Q)-^(Q) 



(16) 



Here La_^(Q) = La — (L^ • Q)Q is the component of L^ 
perpendicular to Q and Fa{Q) is the form factor. 

It is important to discuss how the directions of the 
Pauli-matrices a are fixed in the usual situation in which 
the directions of L are known with respect to the crys- 
talline axes and the difference in the present case. The 
quantization axis of the neutron spin is fixed externally 
to the sample by applying a (small) magnetic field in a 
specific direction, with respect, say to the momentum 
transfer direction Q of the neutron. This fixes a^ with 
respect to the crystalline axes and the experiment is done 
with various choices of 3 with respect to Q. The other 
directions 1 and 2 are then fixed through knowing the 
direction of L with respect to the crystalline directions 
and the use of the dot-product in ([T5|). In effect, B can 
be treated classically in such situations. 

This is to be contrasted with the present situation in 
which the basis vectors of the ground state (±1, ±1) spec- 
ify only the direction of the orbital moment, up or down 
(or zero) as being along the normal to the cu-o planes 



denoted here by the z — axis. Taking matrix elements 
of B in the ground state ^ leads to off-diagonal terms 
in these basis vectors. As noted this cannot be specified 
as a magnetic field operator generated by magnetic mo- 
ments in specific directions with respect to the crystalline 
axes; all that can be said is that the off-diagonal matrix 
elements are matrix elements of a magnetization oper- 
ator orthogonal to the direction z with which the basis 
vector are specified. A purely real ground state wave- 
function means that only matrix elements of S^ and T^ 
generate the off-diagonal elements and so only L^(Q) 
enters in in (1161) . Correspondingly, only cti appears in 
([T6|. There is no way to fix 1 with respect to the crys- 
talline axes. If however one used a more general choice 
of the wave-function so that it is complex, Lx(Q) are de- 
termined by matrix elements of S^,T^ as well as S'^,T'^, 
and correspondingly 172 enters in the calculation. The fi- 
nal answer for the spin-flip cross-sections of the neutrons 
cannot (and does not) depend on the choice of the wave- 
function, nor can the directions 1 and 2 be determined 
with respect to the crystalline axes. 

The difference in our case from the traditional case 
arises from the fact that L does not come physically from 
an atomic orbital moment where the three different com- 
ponents of the orbital angular momentum can be defined 
with respect to the crystalline axes. Rather, in our case 
only the z-component is defined in the basis; the mixing 
in the ground state of the basis is due to a transverse 
field operator as discussed above. The physical basis for 
the transverse field will be specified below. 

Our purpose is to interpret experiments which deduce 
everything from measuring (functions of) the matrix el- 
ements of (t(Q) through three different choices of the 
quantization axis 3 with respect to the Q. From (|16p. 
it follows that when the polarization of the neutron, i.e 
the direction 3 of cr is chosen parallel to Q, there is only 
spin-flip scattering while for any other choice there is 
both spin-flip and spin-nonfiip scattering. We have iden- 
tified L^ as proportional to Lz- As explained, we have a 
freedom of choice of rotating the 1 and 2 directions of the 
neutron a by any arbitrary angle about its chosen 3-axis. 
It thus follows that for any choice of Q and the neutron 
polarization, one can never determine h^ and Ly. One 
can only determine L^ and the component of L perpen- 
dicular to it, which we will call L^. Having shown this, 
we can do the calculation in the simplest choice in which 
the wave- function is real and only the operators S^,T^ 
and K^^ appear in the Hamiltonian. 

The situation may be contrasted with the case when 
the direction of order of the system is fixed by, for exam- 
ple, crystalline anisotropy as in the anisotropic Heisen- 
berg model. The order parameter (M) is then fixed with 
respect to the crystalline axes and may be regarded as a 
classical source for a classical B(r) in Eq. (IT71) . The mag- 
nitudes of two of the components of {'M.x,M.y,'M.z) can 
then be determined by measuring the neutron scattering 
cross-sections by polarizing the neutron beam in two dif- 
ferent directions with respect to a momentum transfer Q 



and using Eq. (TTGI) . One can then change Q and repeat 
the measurement to determine all the three directions of 
(M) for simple magnetic order (or measure at other Q 
for more complex order.) If one has an Ising model in a 
external transverse field, the direction of M is similarly 
fixed. 

The problem discussed above is different also from 
problems of scattering neutrons in the quantum Heisen- 
berg antiferromagnets, where the kinetic energy terms in 
the Hamiltonian are quadratic operators in the spins, as 
opposed to the present case where they are linear. In 
such cases, quantum-mechanics only induces a reduction 
in the ordered spin-moment without changing its orien- 
tation through a Debye- Waller factor due to zero point 
spin fiuctuations while transferring weight to an incoher- 
ent background. 



IV. MATRIX ELEMENTS FOR NEUTRON 
SCATTERING 

We shall show here that there are two kinds of matrix 
elements in scattering of neutrons due to the linear com- 
bination in the ground state ([5]). This can be seen most 
clearly from the rotational and time-reversal properties 
of the basis states written in the form of Eqs. ([6|). The 
basis states are connected though what we might call the 
time-reversal part of the dipole Hamiltonian which in- 
volve L'^a" + L^cr+. This is the usual scattering. But 
as already mentioned the basis states also go to each 
other under successive 7r/2 and 7r-rotations in real space 
through the axis normal to the plane at the center of a 
cell. We show here that the dipole interaction (ITTI) has 
a finite projection to such rotation operators also. We 
will call these matrix elements of the "rotational" kind. 
The matrix elements from a given initial basis state and 
a given final basis state for both kinds must be summed 
and then squared to get the scattering cross-section. 

Since the local loop current is actually an extended ob- 
ject, the neighboring loop currents although orthogonal 
have a finite matrix element through the spatial depen- 
dence of the dipole interactions Hamiltonian: 



Hint (R) 



cr-L-3(cr.(r-R)(L- (r - R) 
|r-R|3 

1 

cr • Vr X [L X Vr 



Rl 



(17) 



Here (r — R) is the unit vector along (r — R). r and R 
is the position vector of neutrons and local moments. 

Let \ipa) denote one of the four loop current state 
I ± 1, ±1) and \ipa, R-i) denote the local moment state of 
ipa located at R^. This local moment state can be written 
as the direct product of the coordinate part and a " mo- 
ment" part as |-0a,R-i) = |0(R- — R-i))IV'o, R-i)s- Here the 
moment part \ipa, Ri)s are the four possible states | ± 1) 
and |0^'^). The coordinate part (/>(R— Ri) describes the 
finite size distribution of the magnetic moment centering 



around R^. For orbital moments, (j) ma-y be taken to be 
the real part of the wave-function while ^ rnay be taken 
to be a phase varying around the loop for the loop cur- 
rents. For 0(R — Ri) = 5(R — Rj), there are only the 
usual matrix elements of the dipole interaction between 
the basis states of [2|) because 

/"d3i?0(R-R,)ff,„t(R)0(R-Rj) 

= /"d3i?5(R-R,)-ff™t(R)(5(R-R,) oc5,,.(18) 

We will call the matrix elements of the dipole interactions 
for (i — j) the matrix elements of the "spin" kind. 

But for finite size 0's, the above matrix element is in 
general not zero. The detailed form of this function 0(R) 
is not very important to us. Here we also assume that the 
distribution function satisfies the orthogonal relations as 
/ d^R(j){Ii ~ Rj)0(R - Rj) = (5.y. 

Then such matrix element are 

(^a,R*|^,„t(R)IV'b,Rj) 
= - f d^R(l>{R-R,)cr-B(f>{R-Rj) 

1 

B = Vr X [Lij X Vr 



Here t„ is the reciprocal lattice vector of the lattice. Sup- 
pose the unit vector of transfer momentum is q, then the 
scattering intensity is 



/(Q)(x|(L^(Q))|^ 



(25) 



Ir-R| 



Here Lij = {ipa,Ri\sMi-'b,'R-j) s- Wc will call them matrix 
elements of the "rotational kind" . 



A. Matrix element of the "spin" kind 

The ground state expectations of L.^ using Eq. ([7]) is, 

(Li(Ri)) ^ (sin6l,cos6i), (19) 

(L,(R2)) -(sin0,O) (20) 

{U{R3))^{-sm9,~cos9), (21) 

(Li(R4))-(-sin0,O), (22) 

where now the three components refer to the "directions" 
t and z, respectively. 

The calculation in Eq. (|19p is only for one unit cell. 
The magnetization in the lattice just repeats the same 
result in each unit cell. 

L(Q)=^L(R,,)e-^Q-^- 

and the magnetization is given by 

Lt{Q) = ^sine'(e-*<^-^i +e-*<^-^^ 

n 

_g-«Q.R3 _ e-*'^-^^)5(Q - x„) (23) 

L,(Q) = J2 cos^(e-'Q-^^ - e-'Q-^^)<5(Q - t„) 

n 

(24) 



Since we treat the orbital current loop as a point-like 
spin, the scattering amplitude is a constant as function 
of Q. So, this calculation is not designed to give the 
structure factor correctly. 

Suppose this is the only contribution to scattering. 
Then at any Q = (a*, 0,0), barring multiple cu-o lay- 
ers per unit-cell (which we will have to consider for 
YBa2Cu20e+s, the "tilt" angle would be determined by 

\Lt\ _ 2 sin 6* • 2sin(Qro) 
\Lz\ cos 6 ■ 2 sin{Qro) 

For the fitted parameters from calculation of the collec- 
tive modes in I, we have sin6' — —0.26, so the tilt angle 
is only about 28°. 



B. Matrix elements of "Rotational" kind 

We now consider the second kind of scattering starting 
from Eq. 1^^. 

We can make use of the fact - = f d^ge*'*'"^ to trans- 
form the above result to momentum space. 

(Va,Rz|ffmt(Rl)|V^b,Rj) 

= cr • /" d^Qe'^'' Q X {Lij X Q)S'*^(Q) (26) 

with S^jiQ) = f d^R(l){R-R,)(t,{R-Rj)e-'^-'^ 

Now we can expand the exponential factor and the lead- 
ing nonzero term is 

StA<i) ~ -i I d^R<i)[^ - Ri)(l){R - Rj)(Q • R) (27) 

To be specific, we consider 6*12 (q) first. For simplicity, we 
assume that the distribution is isotropic (f){R) — (t>{\R\). 
Since Ri = (rg, rg) and R2 = (— tq, rg), then it is easy 
to see that (/)(|R— Ri|)0(|R — R2I) is even function about 
Rx and R^. Thus we have 

5i2(Q) « -iqy f d^Rcb{R-Ri)cb{R-R2)Ry 

= -iQyC (28) 

S23{Q) ~ -iqx fd''Rcf>{R-R2)HR-R3)Rx 

= iQrC (29) 

534(Q) « ^iq, f d^Rq){R-R3)(l){R-R4)Ry 

= iQyC (30) 

SiiiQ) ~ -iqx f d'^R(biR-R4)(j)(,R-Ri)Rr 

= -iQ.C (31) 



where C = fd^R(j){\Il - Ri|)0(|R - RaDi?^- We see 
that C is of O(ro). We will find the coefficient only by 
fitting to the data. 

Note that this new term is also proportional to Lx(Q), 
just as the traditional matrix elements for polarized neu- 
tron scattering. But the momentum dependence is dif- 
ferent from the classical neutron diffraction expression. 
We can repeat this calculation at all possible neighbor- 
ing local moments and find similar expressions. Recall 
that the 4 states are collections of 4 local spin 1 states. 
The ground state is superposition of the 4 states, thus it 
is also a collections of 4 local spin 1 states. To compute 
the total neutron scattering amplitude, we need to evalu- 
ate the matrix elements of L between all the neighboring 
local spin 1 states of the ground state. The 4 local spin 1 
states of the ground state is given by (Here we only write 
the spin part) 

|G,Ri), = c2|l)+sc|0i)+sc|02)+sVl) 

|G,R2). = sc\l)+c^\0^)+s^\0^) + sc\-l) 

\G,R3)s = s2|l)-sc|0i)-sc|02)+c2|-l) 

|G,R4),, = sc|l)-s2|0i)-c2|02)+sc|-l) 



Here c = cos9;s — sm9. Then it is straightforward to 
find the following matrix elements for L 



/(f + sin^e) i ^ sin(26l)\ 
) , = , — cos y, 

' V 4 2 4 y 



(G,Ri|,L|G,R2 



(G, R2I.LIG, R3)s = (^ cos2 9, -^ cos9,~^sm{29) 



(G,R3|sL|G, R4 



(1 + sin 



i sm{29)\ 

, - cost', 

4 '2 ' 4 / 

(G, R4|sL|G, Ri). = ( - ^ cos^ 9,-'- cos 9, ^ sm{29) 

Other matrix elements are the complex conjugate of the 
above equations. If we put all the terms together, one 
can see that the imaginary part cancel out and only the 
real part contributes. We find the scattering amplitude 
in the momentum space as 

{G\Hm(R)\G) = 2^cr • Q X iUj X Q)5,j(Q) (32) 
(y) 

with (a,fe) = (1,2), (2, 3), (3, 4), (4,1) and L,, = 
Re(G,Ri|sL|G,Rj)s. Putting all the above results to- 
gether, we find 

(G|i/,„,(R)|G) = -AtQyCa ■ Q x (L12 x Q) 

-AiQ^Ca- ■ Q X (L41 X Q) (33) 

Again we can consider the transfer momentum Q — 
(a*, 0,0). The scattering amplitude is 

(G|iI,„t(R)|G) = -AiQ^Ca- ■ Q x (L41 x Q) (34) 

which can be rewritten as 

(G|iJ,„t(R)|G)=^-Qx[LeffxQ] (35) 

Lcff = 4^Q,G(icos2 0,-isin(20)) (36) 



C. Total Matrix Element 

The form, Eq. psp . has the same structure as the ma- 
trix element of the " spin- kind" Eq. (f25|) . We can combine 
them together to find the total magnetization as 



Ltot — Lo + Lcff 

Lq = 2isin(Qj,ro)(2sin6',cos( 



(37) 



Lcff' ^4iQ^C{-cos^ 



1 



sin(26l) 



where the two components again refer to t and z respec- 
tively. Here Lo is the magnetization form the " spin- type" 
contribution. 

The parameter tq specifies the location of the mo- 
ments. For the ground state determined from the col- 
lective modes, we have sm9 — —0.26. If we assume that 
Q^C ~ —0.5 sin((5a;ro), then we find that for transfer mo- 
mentum Q = (a*, 0, 0) the tilted angle is about 50°. The 
experimental results are 55 ±7° and 35 ±7°, respectivelj^ 
for YBaCuO%,Q. But most of the data available for the 
collective modes from which the angle 9 is deduced is in 
Hgl2Ql with a. Tc ~ Q\K^ which from the general phe- 
nomenology has properties close to those in YBaCuOe,^. 
Collective modes have also been found in YBaCuOQ,gr^ 
at similar frequencies but detailed information about the 
dispersion of the two branches of collective modes is un- 
available. The tilt in Hgl2Ql at (1,0,1) is deduced to 
be^ 45 ± 20 within the estimates provided here. We can 
claim that with G = O(ro), we get the correct trend and 
the magnitudes of the tilt in agreement with experiment 
within the stated error bars. It should be stated that a 
simplified model is used to calculate the collective modes 
and the ground state from which 9 is determined. It 
could easily be ±10° from that deduced. 



V. CALCULATION INCLUDING FORM 
FACTORS 

A. Real space representation of the Operators 

To calculate the neutron scattering including form fac- 
tors, one needs the wave- function which have the orbital 
moments and the microscopic representation of the op- 
erators which lead to the produce the admixed ground 
state wave-function ^. The approximate space repre- 
sentation of the basis states is relatively straight-forward. 
One may represent them with complex wave-functions 
carrying currents around the indicated loops as in Fig. 
(HI). The subtle issue is in the representation of the op- 
erators {Sf + Tl), which give the linear combination of 
the classical basis states in the ground state in ^. They 
were introduced at a formal level in I simply because 
such terms are allowed and because they are necessary 
to calculate the observed collective modes in the loop 
ordered state. But what is the physical origin of such 




Figure 2: Current pattern of the rotation operator in Eq (|38l 



quantum terms? The physical origin was already derived 
in connection with calculating the spectra of the collec- 
tive fluctuations in the quantum-critical regime of the 
loop ordered phase^^. A simpler derivation is given in 
the Appendix |^ The operator Sf + T} is given by the 
current operator schematically shown in Fig. ([2]) and 
its hermitian conjugate which has current flowing in the 
opposite direction. In the continuum limit, this is sim- 
ply the operator proportional to "^^^n^ V x j(r), with r 
measured from the center of each unit-cell and with cut- 
off at the boundary. We recall from I that the operator 
{S} + Tl) acting on any of the four states (±1, ±1) in a 
unit-cell admixes the state rotated by ±7r/2 to it. This 
operator corresponds to the collective part of the follow- 
ing fermion operators (See Appendix VK\ 



E' 



,ikRi 



2i 



. SyC; 



y^x 



SxC' 



^tlLtJ 



^0 



Wx,kPy,k-PyMP^,k) (38) 



Here E^ is a normalization of the wave-function and Px^ 
and Py k are the fourier transform of the oxygen p-orbital 
operators in the x-direction and y-direction, respectively 
around the Cu-atom in the unit-cell i. Such an operator 
is created microscopically through expressing the nearest 
neighbor interactions in each unit-cell in terms of cur- 
rent operators and combining such current operators to 
form closed loops within a unit-cell. Five different point 
group symmetries result including that of ([55]) sketched 
in Fig. ([2]). That this operator admixes the basis states is 
shown in detail in Appendix [X] 

Given this, we can calculate the form factor of neutron 
scattering by further approximations to represent the col- 
lective states shown in Figures ([J) and ©.These are de- 
scribed in Appendix [C] We now proceed with the com- 
parison with the experiments in evaluation of the neutron 
scattering using these results. 




k-Jc' 



Figure 3: M^ (black curve) and Mz (red curve) as a function 
of fez. 



B. Results for tilt and form factor 

It easy to see that if we include the form factor ef- 
fects, Lz decreases faster than Li as Qz increases. (See 
Fig ([3])). (Recall that the subscript 1 in Li stands for 
whatever component of L perpendicular to the z-axis is 
being measured.) We therefore find the tilt angle for 
Q = (101) is larger than Q = (100), qualitatively consis- 
tent with the experiments^. How the tilt angle depends 
on the momentum is determined in our calculation by the 
choice of parameters such as r^ and Tq, which is hardly 
definitive. 

For example, as a reasonable guess, we can take tq — 
0.25a, r„, = 0.175a. Here a is the lattice constant of xy 
plane. The tilted angle is quite sensitive to the current 
width parameter w. From experimental data, we know 
that the intensity of neutron scattering of L = 2 is half 
of that for L = 0. For a Gaussian shape dependence of 
Qz, we can deduce that w « 0.5a. The half width of the 
current is comparable to the radius of the current. 

For the above choice of parameters, the intersection of 
these two circles are Ria = (0.238a, 0.075a) and Rk, — 
(0.075a, 0.238a). If we use the approximate form factor of 
Eq. (|G9p . then for transfer momentum Q = ^(0,1,0), 

we find that the tiled angle is </> = arctan(-U^) « 64°. 

Since the width w = 0.5a is not small compared to a, 
we can not use the above approximation. Therefore we 
have to numerically transform Eq. (|C6|) and (|C4|) to mo- 
mentum space. For Q = ^(0, 1, 0), we find that the tiled 
angle is ~ 26°. If we also include the "rotation" kind 
contribution to the moment from Eq. (|C10|) and (|C11|) 
and choose the parameter C = —0.2a, we find the total 
tiled angle is </> r; 38.5°. For YBCO, c « 3a. Therefore 
for transfer momentum Q = (0, — , — ), we find the tiled 
angle (/) « 47°. 

For YBCO, there are two copper-oxygen planes. Sup- 
pose the distance between these two layers are d, then we 
will have an extra structure factor cos{Qz^)- This makes 
the amplitude oscillate with Qz- 



VI. CONCLUDING REMARKS 
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In this paper we have used the result derived earher in 
I that the ground state of the loop order state is a linear 
combination of the classical basis states. This necessi- 
tates a new view of calculating the neutron-scattering 
cross-section, in which one must consider the flip in the 
neutron spin due also to the matrix elements of flip oper- 
ators of the quantum-moments in the ground state. The 
simplest (mean-field) ground state was considered. This 
does not include the zero-point deviations due to spin- 
waves. We have calculated these to be a small effect 
because the collective modes are all at finite energy. 

Our method of calculating the neutron scattering 
cross-section may be useful in other quantum problems, 
for example the transverse field Ising model, provided the 
transverse field is not an external field with a specified 
direction but internally generated. It may also be useful 
when scattering experiments are done in other loop or- 
dered states such as the anomalous Hall effect and the 
topological insulator states. 

We have been able to show reasonable consistency of 
the calculation with the experimental results but a pa- 
rameter C was introduced which we find hard to calculate 
but can argue only that it should be small compared to 1, 
as is indeed found. Especially gratifying is that the title 
angle deduced depends on the Bragg-vector. This would 
not happen in the traditional usage in which the neutron 
spin fiips quantum-mechanically of a spatially dependent 
classical magnetic field due to the magnetic order. 

It is worthwhile commenting on earlier attempts to 
understand the "tilt". One was based on spin-orbit 
scatteringii. This had two difficulties; there is no such 
term in Hgl201 and the magnitudes do not come out 
reasonably without giving a scattering at (2,0,0), which 
is not in agreement with elastic experiments. The other 
idea is that there is a moment on the triangles made 
through apical oxygens'^^. This is allowed by symmetry®. 
Mean-field calculations do not provide any significant 
such moment for any reasonable set of parameters. There 
is a more basic problemm with having any significant am- 
plitude for such a contribution. It is that such moments 
would provide zero contribution of the moment perpen- 
dicular to the z-axis for any scattering vector with Q2 = 
due to the structure factors in the cu-o lattices with a 
plane of reflection in the unit-cell perpendicular to the 
z=axis. A large tilt has been deduced for Q = (1, 0, 0). 

Our calculations in this paper also have an impact for 
some other measurements. If the tilt were due to mo- 
ments in triangles made through apical oxygens as dis- 
cussed in the previous paragraph, there should be mea- 
surable magnetic fields— at several sites in the lattice, 
detectable in NMR experiments. With the quantum ori- 
gin of "tilt fields, the calculations based on a classical 
magnetic field tilted in specific directions are not valid. 
Further consequences of the quantum origin of the "tilt, 
for example in NMR experiments is well worth further 
investigation. 
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Appendix A: Microscopic Theory of the kinetic 

energy operator (S^ +T^ + S'^ + T^) for the Loop 

current states 

The loop current states were derived^ starting from 
the basis of the three-orbital {di,pix,Piy) per unit-cell. 
Besides the kinetic energy operator between the neigh- 
boring d,px and py orbitals, the model includes local 
repulsion at each site, (the U's) and near neighbor re- 
pulsions of charges: VniUj. The on-site repulsions are 
assumed to serve only to renormalize the remaining terms 
in the Hamiltonian in the metallic state. A crucial role is 
played by the operator for nearest neighbor interactions. 
For spin-diagonal terms in the operator riirij, it may be 
written as 



^E- 



t,cr"'j,<T 



= V/2{\j,,\^-n,-n,), 



(Al) 



Here jij = i ^^ c^a^i-'^ + ^•'^- ^^^ ^^^^ operator content of 
a current. Discarding the one particle terms, a mean-field 
approximation is made on |jijp: 



\Jv\ 



- rjtj 



* jij + fiuctuation operators) (A2) 



Here rij is the collective part of jij . Now the effective ki- 
netic energy on the link (ij) is tij +iVrij /2 which gives a 
complex kinetic energy with a phase « rVij /{2tij). Phase 
differences on links within a unit-cell are combined to 
form closed loops which have invariant fiuxes with differ- 
ent point group symmetries. For non- intersecting loops 
on the cu-o lattice, there are five and only five such closed 
loops possible^-- . Two of these which transform as dou- 
bly degenerate vectors Ex , in the {x ± y) directions and 
their time reversed partners form the four flux patterns 
depicted in Fig.([T]). Simple mean- field calculations as 
well as more elaborate calculations^^ show that this is 
the most stable allowed symmetry and experiments^ have 
shown it to be the realized state in all the under-doped 
cuprates investigated. 

The other current loop which plays a crucial role in our 
considerations here and in the collective modes of I and 
played an important role in the calculation of collective 
fluctuations^® in the quantum-critical regime is the pat- 
tern transforming with the full symmetry of the lattice. 
It is depicted in Fig. ([2]). We will show that it has local 
matrix elements between the states (±1,±1). Before we 
show this it is important to specify some properties of 
the basis states (±1,±1). 
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Corresponding to each of the four collective states 
(±1,±1), there are four kinetic energy Hamiltonians for 
the fermions: 



iJ(k) = 






'^tpd^x 


'^tpd^y 


'^f^pd^x 





^pp^xSy 


ItpdOy 


TppS'xSy 






(A3) 



with 5a; = sin(A:a;a/2 + (j)x), Sx = sin(fcj,a/2 + (f>y), 
Sx = sin(fca;a/2), Sy = sin(/cya/2). Then the four dif- 
ferent mean field Hamiltonian correspond to the order 
parameter {(j)x,4>y) — (±J7,±J7), where fi is the magni- 
tude of the order parameter which is determined vari- 
ationally. The eigenvalues and eigenstates of the three 
fermion bands in the zeroth order of tpp/tpd are given by 



Ey, = tpd\J SI -I- 52 

V 2 ^ Oxy Jxy ^ 

|2,fc) = (0,-|^,|^)^ 

^ ^xy ^xy ^ 

|3,fc) = — fi,— ,^)^ 

V 2 ^ ^xy ^xy ' 



Sx + s,,. 



with Sxy -- y^x ' -y 

To find the collective operator which rotates among the 
four states (±1,±1), let us first find the operator which 
rotates among the fermion states of these four collective 
states. Then the collective operator we are looking for 
is simply the collective (but uncondensed) part of such 
an operator formed from the fermions. Let us use sub- 
script a, b, c, d to label order parameter (1, 1), (—1, —1), 
(—1,1), (1,-1). Then we can find a unitary operator 
which rotates among the four states. 

Uad = |l)d(l|a + |2)d(2U + |3)d(3U 

1 



with Sx 




1 



-2nsxCy/sly 1 



2nsxCy/sly +o{L) 



X ' "y 



Sy. Similarly, we have 



Uac = |l)c(l|a + |2)c(2|a + |3)c(3|a 

1 

= 10 1 



2Q,SyCxlsl 



ZllSyCx/ S,j.y 

1 



,{L) 



We can now introduce the generator of the above trans- 
formation 

aik = '^^^iplkPv.k " plkP^,k), (A4) 



E, 







4,k = -^^^^ipl,kPy,k - pl,kPx,k) (A5) 



In terms of the generators, we have 



Uac 



1 



^^ai. 



fc' 



Ucb 



1 



ing^ 



Ubd = I - i^g: 



/ 



,fc' 



f/da = 1 - i^g^ 



y,k^ 
y.k 



/ 



J 



It follows that the fermion Hamiltonian, ^j^ 9x k ~^ 9v k 
serves to mix any of the fermion states corresponding 
to the collective states (±1,±1) with a state rotating ft 
by ±7r/2. The current corresponding to this term runs 
within the unit-cell from one oxygen to the next in the 
clockwise direction as shown in Fig. ^ and the corre- 
sponding current in the anti-clockwise direction. It has 
already been shownlS that such a term is generated by 
the nearest neighbor interactions. It then follows that a 
collective current state of the same symmetry also exists 
which serves as a kinetic energy term mixing the collec- 
tive configurations. In the basis of the collective states 
I ± 1, ±1), this has all the transformation properties and 
commutation rules of the operator {S^ + T^ + S"^ + T^). 



Appendix B: Need for a 4-dimensional 
representation of the operator S 

As discussed the gauge transformation is given by the 
rotation matrix 



U{e) = e*''i'S''/2g*e2TV2g»e3-R-^V2 



(Bl) 



sin ( 



U{e)K^^U{ey =cos6iicos6i2/s: 



The classical AT model is invariant under this transfor- 
mation. Under this rotation, the quantum term S^, T^ 
and K^^ will be transformed as follows 

U{e)S^U{e)^ = cos 6(3(008 6*1 51 - sin6iiS'2) 

-sin6l3(sin6'iXi3-hcos6liis:23) (b2) 

U{9)T^U{9y = cos6l3(cos6l2Ti - sineiaT^) 

6'3(sin6'2X^^ +cos6l2X^^) (B3) 

^^o\t ^ COS 01 cos 6I2/S:" + sin 61 sin 6*2 if ^^ 
(cos 01 sin e2K^'^ + sin Oi cos 6*2^4:^^) (B4) 

The quantum terms related by this transformation are 
equivalent to each other. In the mean time, the mean 
field ground state is also transformed to a new form. In 
the loop current state basis, the original ground state can 
be written as \G) — {c^ , sc, sc, s'^)'^ . After the transfor- 
mation, we have \G') — U{9)\G) which is a superposition 
of the 4 states with complex number coefficients. Thus 
after the transformation, (S^) is nonzero and depend on 
the gauge parameter ^1,2,3- 

Since M(q) as a physical observable is gauge invari- 
ant, it should not depend on the gauge parameter 0i,2,3- 
This gauge dependence should disappear if we make a 
corresponding unitary transformation on the local spin 
operator S^''^'^. But now we meet a immediate difficulty. 
The transformation U{9) is a 4 by 4 matrix simply given 
by 



C/(6')= diag{e'*'^e"^^e 



'} 



(B5) 
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with 



1+92+93 



1+^2 



and above function. If the current width is narrow, we have 



91-92+^3 



Translating this transformation to the 
basis of local spin 1 state, we find that the corresponding 
transformation is 



V{d) = diag{. 



,i't>l „i4>2 



+ e 



i<p3 „i<t>i 



} 



(B6) 



Clearly, on can see that V{6) is not a unitary transforma- 
tion. Therefore we cannot use V{9) to cancel the gauge 
dependent. 

This difficulty arise because there are 4 states in the 
loop current state basis but there is only 3 states in the 
local spin 1 state basis. In order to overcome this diffi- 
culty, we have introduced a 4 dimensional representation 
of the spin 1 states given in the paper. 



Appendix C: Approximate Calculations of Form 
factor 



We first consider the form factor for the Lz- Since 
Lz only has diagonal matrix elements, we do not need 
to use the rotation operator yet. The form factor is de- 
termined by the current density around the 0-Cu-O tri- 
angle. To simplify the calculation, we approximate the 
triangle shape by a circle shape current density and also 
assume the density is a Gaussian distribution around the 
circle. Then, the wave-function of the local spin states 
in terms of cylindrical coordinates r = (r cos (j), r sin 0, z) 
are given by 



V'a, 



1 



±1 



-,±i4 



2ttw 



V'a,0 = 



2ttw 



(CI) 
(C2) 



Here we assume the radius of the circle is r^ and cen- 
tered around Ra and w specify the width of these circles. 
Now we can verify that ^pa.i gives a torus shape current 
density. 



j(r) 



2ttw 



2 



Z72 — 



(C3) 



Similarly, ip^_ gives a torus shape current density with 
opposite direction. -00 is real function, so the current is 
zero. These 3 states correspond to | ± 1) and |0) point 
like angular momentum states. If we take w ^>- oo limit, 
we find a current circle j(r) = S{r — ro)S{z)(p. 

We still have L^ = ijj'^Lz'4', where L^ = —i-ffz- The 
(j) dependent part will give the same matrix elements as 
before. Now we also have an extra r and z dependent 
part as follows 



Here r„ = ^/(a^ 



1 



2ttw 



■ exp 



-2wl 



2w2 



Ra,xY + {y- Ra.yY with fl 

Then the form factor is the Fourier transformation of the 



(C4) 
1,2,3,4. 



F(k) = / F(r)e -'•'•'■ dV « f{kr)e' 



ik R„ 



with 



cell 

fiK 



1 



y2TTw Jo 
Using this result in Eq. ([24| . we find 



' Jo{krr)rdr 



i.(k) 



COS! 



F,a(k)e-^'^-^i-F,,i(k)e- 



ik-Ra 



(C5) 



The form factor of Lt is more complicated, since the 
rotation operator has its own coordinate dependence. As 
explained before, the rotation operator is currents fiows 
around the four oxygens in the unit cell. Here for sim- 
plicity, we approximate the current around the oxygens 
as a big circle with radius ro centered at R°. We have 
Lt = il^'^Lt'ip- Here Lt will annihilate or create the e"'^ 
factor. The (p dependent part will give the same matrix 
elements as before. For the local spin state at Rq, we 
also have an extra r and z dependent part as follows 



FxA^) 



1 



2'KW'^ 

X exp 



exp 



1 z^ 

2w^ 2w'^ 



1 

2w 



^{ra 



rsf- 



z 
2w^ 



(C6) 



Here r — \/ x^ + y^ . It is easy to transform the z de- 
pendent part to momentum space. The a:, y dependent 
part has to be computed by numerics. Including the form 
factors in Eq. (j23l) , the x-component of moment is 



i*(k) 



sm( 



F,,i(k)--*-^i 



-^x,3(k)e 



-ik-Ra 



F.,2(k)e-' 
F,,4(k)e-*-^' 



kRa 



(C7) 



If the current width is quite narrow w ^ a, here a is 
lattice constant of xy-plane, then Eq. (IC6p can be ap- 
proximated by to two Gaussian peaks located at R;!^ and 
Rq where R^ and R^ are the two intersection points of 
the two circles. Thus we have 



^x,a(r) 



2'Kvr- 



(C8) 



Transferring to momentum space, we find the form factor 



^a;,a(k) = yfUn 



■ (e"'-^" 



ik-Ri 



(C9) 



with a = 1,2,3,4. 

So far we only considered the "spin" kind contribution 
to the moments. Similarly, we can also compute the form 
factor for the "rotation" kind contribution. Follow the 
same line of arguments, the (j) dependent part will give 
the same matrix elements as before. For the the matrix 
element {ipa\Lt\tJJb), there is an extra r and z dependent 
part as follows 



12 



^x,ab(r) 



1 



X exp 



exp 



^^'-''^'-ii;^ 



exp 






2 



1 ^ n2 ^' 



z2 - 

Similarly, or the the matrix element (i/iajL^I?/;;,), there is an extra r and z dependent part as follows 



F,^ah{r) 



1 



2t^w- 



■ exp 



1 z^ 

——{ra-r,f- -— 



exp 



^11)^ Aon^ 



2 

4w2 



(CIO) 



(Cll) 



And the form factor is the above function transform to momentum space. Therefore, for transfer momentum k = 
— (0, 1,0), we find the effective moment due to the "rotation" kind contribution is 



-eff = -«fcy^(^-,i2(k)(l + sin2 0), 0, F,s2{^)sin2e 



(C12) 



These results have been used in the paper to calculate the form-factor and the tilt angles and compared with 
experiments. 
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